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1.  BACKGROUND  AND  PURPOSE 
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It  has  been  found 11  that  the  probability.  Pd.  of  achieving  a  given  degree  of  damage  to  a 
point  target  is  a  log-normal  function  of  the  environment  intensity.  E  (total  dose,  in  rads  tissue, 
for  example),  produced  by  the  nuclear  burst: 


where  erf(.v)  is  the  error  function  and  E-M  and  <r  are  the  two  parameters  that  characterize  the 
vulnerability  of  the  given  system  to  the  environment  in  question.  The  vulnerability  parameter  £'50 
is  that  value  of  the  environment  intensity  that  produces  a  value  of  P„  =  0.5.  and  a  is  a  measure 
of  the  steepness  of  the  P„{E)  curve.  Thus.  <r  and  Ew  depend  on  the  physical  nature  of  the  system 
as  well  as  on  the  identity  of  the  environment  that  induces  the  damage — neutron  fluence;  peak 
gamma  dose  rate;  total  dose;  blast  peak  static  overpressure,  impulse,  etc;  and  electromagnetic 
pulse  (EMP). 

Present  methods  for  treating  extended  targets  require  either  point-by-point  calculation  of 
damage  probability  for  each  point  target  within  the  extended  target  (see  below)  or  that  the  analyst 
resort  to  existing  manuals1-'  that  deal  with  specific  strategic  targets  and  tactical  troop  and 
equipment  responses  to  the  burst. 

For  new  extended  targets,  or  for  extended  targets  that  have  not  been  previously  treated  in 
these  manuals,  the  calculation  of  the  damage  to  an  extended  target  (composed  of  a  number  of 
discrete  point  targets  and  of  finite  extent)  that  is  exposed  to  nuclear  environments  requires  a 
weighted  average  of  the  P „  values  representing  each  discrete  point  target.  This  calculation  can 
prove  excessively  lengthy  for  those  extended  targets  that  are  composed  of  many  point  targets; 
indeed,  assigning  appropriate  weights  is  a  problem  itself.  These  calculations  are  lengthy  even  if 
one  resorts  to  nomograms  or  a  slide  rule  that  have  been  developed  to  speed  the  calculational 
process.1*-7  It  would  be  convenient  to  be  able  to  treat  such  an  extended  target  as  a  single  target 
that  requires  a  single  P „  calculation. 

The  purpose  of  this  study  is  to  develop  techniques  that  mathematically  treat  extended  targets 
as  point  targets.  The  analysis  develops  appropriately  accurate  expressions  for  calculating  the 
damage  to  extended  targets  and  methods  for  reducing  the  number  of  required  calculations. 


1  W  I..  Vault  an<J  W.  I-..  Sweeney.  Vulnerability  Data  Array  Progress  Report  KY7h.  PY7T  (Cl.  Harry  Diamond 

1  aboratories.  HDI.-PR-77-t  (December  IV77|.  (CONI  ) 

‘  U  1  Dtirgin.  D  R  Alexander,  and  R  N.  Randall.  Hardening  Options  for  Neutron  effects.  f  inal  Report.  BDM 
Corporation.  Harry  Diamond  Laboratories.  HDL-CR-74-052-1  (15  November  |V7fo. 

’  K  N.  Stevens.  Neutron  failure  fluence  Distribution  for  Semiconductor  Devices.  II  I. I  Annual  Conference  on 
Nuclear  and  Space  Radiation  effects.  Seattle,  WA  (24  to  27  July  1772) 

1  Department  of  the  Army,  Staff  Officer's  field  Manual.  Nuclear  Weapons  employment  (Cl.  field  Manual  101-41- 

2  (fehruary  I  Wife  (SP.CRLTl 

'Defense  Intelligence  Agency.  Physical  Vulnerability  Handbook — Nuclear  Weapons  (l  i.  AP-SMI- 1 -2-1  NT  (June 
IVhd)  (CONI  I 

*  C  Stuart  Kelley.  Stacey  I  (iehman.  John  If  Wasilik.  and  William  D  Sc  half.  Nuclear  Damage  to  Point  Targets. 
Harry  Diamond  l  aboratories.  HIM .- I  R- lX7h  (December  I77X) 

’  C  Stuart  Kelley.  John  H  Wasilik.  and  William  D.  Scharf.  Supplement  to  Nuclear  Damage  to  Point  Targets  (Cl. 
Harry  Diamond  Laboratories.  HIM-  PR  I S7P-S  (November  177X1  (SM  Rf  Tl 
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2.  METHODS  FOR  CALCULATING  DAMAGE  TO  EXTENDED  TARGETS  COMPOSED  OF 

DISCRETE  POINTS 

To  perform  an  average  of  a  number  of  values  of  Pd  using  equation  ( I ),  the  analyst  adds  terms 
of  the  form  erf(.r , )  +  erf(jr2)  +  .  .  .  .  Unfortunately,  a  basic  calculational  problem  arises  that 
involves  the  presence  of  the  error  function  in  equation  (I).  The  error  function  is  not  a  closed- 
form  expression  and  must  be  evaluated  either  by  reference  to  a  table  or  by  the  use  of  a  series 
expansion."  It  obeys  no  simple  mathematical  properties  of  addition.  Thus,  there  is  no  simple  way 
to  express  erf(jr,)  +  erf(jr2).  an  operation  that  is  needed  for  the  treatment  of  extended  targets. 
The  remainder  of  this  section  deals  with  approximations  that  can  be  applied  to  circumvent  this 
problem. 


2. 1  An  Approximation  of  Pd 

The  main  span  of  interest  to  the  analyst  in  damage  calculations  is  0.1  <  P ,,  <  0.9  because, 
for  the  tactical  analyst,  damage  probability  less  than  0.1  implies  almost  complete  effectiveness 
of  the  unit,  while  Pd  >  0.9  implies  almost  complete  ineffectiveness  of  the  unit  in  carrying  out  its 
mission.  The  FM- 101-3 1  manual  series.*  in  fact,  considers  values  of  just  P ,,  ~>  0.3  as  operationally 
useful.  Likewise,  a  probability  greater  than  0.9  that  an  item  has  suffered  a  given  degree  of  damage 
is  adequate  characterization:  resolution  as  to  whether  P.,  is  0.93  or  0.94  is  not  necessary.  There 
are  occasions,  such  as  in  sensitivity  analyses,  where  altering  an  input  parameter  of  the  scenario 
may  shift  the  P ,,  value  from  inside  to  outside  the  range  of  0.1  to  0.9.  In  such  cases,  the 
approximations  (and  their  limits)  developed  as  follows  should  be  kept  in  mind. 

We  define 


*  *  — i — fn  (E/Em)  .  (2 1 

X  2  <t 

so  that  equation  (1)  becomes 

Pd  =  {[•  +  crf(Jtr)  ]  (3 1 

From  figure  I  it  can  be  seen  that  over  the  range  of  interest  (0.1  <  P ,,  <  0.9.  corresponding  to 
-0.8  <  x  <  0.8).  the  error  function  may  be  replaced  by  a  linear  function.  A  least-squares  fitting 
process  (from  John  Wicklund  of  the  Harry  Diamond  Laboratories)  gives  a  result  that  is  (within 
3  percent)  the  same  as 

erf(.tr)  x  .  (4) 

which  is  depicted  in  the  figure.  Over  the  range  of  interest,  x  approximates  erf(.v)  to  within  a 
fractional  error  of  10  percent.  With  the  approximation  of  equation  (4).  P,t  becomes 

^  +  ~  *-  fn (£/£-,„)  (5) 

2  2v  2  <r 


*  IX-partmcnt  of  the  Army.  Staff  Officer's  field  Manual.  Nuclear  Weapons  f.mploy  mem  ill.  field  Manual  HIM  I 
:  ( February  IWi  (Sl.l'RKT) 

"  M,  Ahramowit/  and  I.  A.  Stcgun.  editors.  Handbook  of  Mathematical  functions  with  formulas.  Graphs,  and 
Mathematical  Tables,  National  Bureau  of  Standards.  Applied  Mathematics  Scries  S' (June  I9WI. 
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erf(x) 


Figure  I.  Linear  approximation  to  error  function:  crf(— at)  =  -erf(.r). 


and  is  accurate  to  within  a  10-percent  fractional  error  over  0.  J  <  Pd  <  0.9.  This  fractional  error 
is  consistent  with  errors  in  the  data  to  which  the  algorithms  below  were  fitted  to  the  environmental 
intensity,  £(r).  Clearly,  we  need  to  carry  calculations  to  no  higher  accuracy  than  that  character¬ 
izing  these  data. 

We  can  characterize  all  the  environment  intensifies  as 

E  =  ArBe<r  (6) 

with  a  high  degree  of  accuracy.®  The  uncertainty  of  the  E(r)  data  to  which  we  fit  these  is  accurate 
only  to  within  a  fractional  error  of  10  percent,  sometimes  only  to  within  30  percent.10  These 


*  W.  K.  Sweeney,  Jr..  C.  Moated,  and  J.  Wickluml.  Nuclear  Weapons  environments  for  Vulnerability  Assessments 
to  Support  Tactical  Nuclear  Warfare  Studies  (U).  Harry  Diamond  Laboratories,  HDI.-TM-77-4  (June  1977).  (CONI') 

10  I..  G  Mooney  and  R.  L.  French.  Improved  Models  for  Predicting  Nuclear  Weapons  Initial  Radiation  environments 
(U).  Radiation  Research  Associates,  Inc.,  DASA-2615.  RRA-T9J  (}|  December  1969)  (SRD/CNWDI) 
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algorithms  were  developed  in  order  to  provide  simple,  yet  accurate,  equations  for  calculating  the 
nuclear  weapon  environments  (static  overpressure,  total  radiation  dose,  thermal  radiation,  neu¬ 
tron  fluence,  etc.).  The  constants  A,  B,  and  C  vary,  of  course,  from  environment  to  environment. 
The  constant  A  depends  on  weapon  yield.  Thus,  the  nuclear  environments  depend  solely  on  the 
weapon  yield  and  the  distance  from  the  burst. 

2.2  The  Effective  Environment  Intensity 

The  approximate  form  of  equation  (5),  has  useful  calculational  consequences  that  will 
now  be  developed.  Consider  an  extended  target  consisting  of  m  different  types  of  point  targets 
(trucks,  people,  radios,  etc.),  each  type  denoted  by  the  subscript  j.  By  type  of  point  target  is 
meant  all  items  that  have  identical  crj  values  and  identical  Ei0j  values.  Let  tt,  denote  the  number 
of  point  targets  of  each  target  type  j,  and  let  the  subscript  i  identify  the  particular  item.  Then  the 
total  number.  A/,  of  point  targets  within  the  extended  target  is 

m 

2  •  (7) 

)=  1 

We  now  calculate  the  value  of  Pd  appropriate  to  the  extended  target  for  a  single  category  of 
damage.  As  the  first  step  in  this  calculation,  we  assume  that  all  point  targets  characterized  by  the 
same  j  subscript  have  the  same  "military  value."  The  /'th  portion  of  the  extended  target  is 
composed  of  n,  items  (each  of  these  items  is  characterized  by  the  subscript  /  and  is  thus  identified 
by  the  subscripts  /  and  /').  We  define  the  P ^  value  of  the  extended  target  as  the  average  of  the 
Pdji  values  for  the  individual  points.  Thus,  P ^  is  given  by 

1 

P lit  ~  ~  2  P <lji  (8) 

ni  (-1 

Therefore,  Pdi  represents  the  probability  of  damage  to  the  /'th  portion  of  the  extended  target. 
Using  equation  (5), 


P,u  *,{  2  2/n<£*/£*w>}  •  W 

where  Eu  is  the  environment  intensity  at  the  /th  item  of  the  /'th  type.  We  replace  Em,  by  E,0J 
because  each  of  the  /th  items  have  the  same  E50;,  value.  Equation  (9)  can  be  arranged  as 


P«  = 


2vb- 


(10) 


If  we  now  define  Erj  to  be  the  "effective  environment  intensity"  of  the  /'th  portion  of  the 
extended  target  by  the  relation 


Pe,  -  ft  £]/"'  •  (ID 

i  =  I 

we  observe  that  the  expression  for  Pm  is 

Pdi  ~  2  2\^2  <T  ^(P-ri/E )  (12) 

Thus,  because  the  approximation  erf(x)  =  x  [equation  (4)J  is  accurate  for  cases  of  interest 
to  the  tactical  nuclear  warfare  analyst,  we  can  avoid  the  necessity  of  using  series  expansions  or 
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tables  to  calculate  erf(x).  This  approximation  led  us  to  the  concept  of  the  effective  environment 
intensity.  This  concept  has  important  and  useful  mathematical  properties  that  are  described 
below. 

2.3  Incorporation  of  Shielding  Factors  and  Military  Values 

We  now  redefine  the  number  of  target  types  by  grouping  those  having  identical  shielding 
factors  and  identical  military  value  in  addition  to  having  identical  vulnerability  parameters  (o- , 
and  £50j).  The  shielding  factor,  Tj,  is  the  ratio  of  the  shielded  environment  intensity  to  the 
unshielded  environment  intensity.  Then,  using  equation  (12).  one  can  write 


Frf/shielded)  =  -/  1  H — j= —  fn  (T 
2  I  V2  er 


jE  rj/  £50j)| 


2  f  '  +  \/2  a  fn{E^/E^  +  ^  ^  <fn(^)| 


which  reduces  to 


£, //shielded)  =  /’^/unshielded)  + 


en(Ts) 


2v2  o-j 

Conversely,  one  can  define  a  shielded  effective  environment  intensity 

E,, /shielded)  =  r;£,,j(unshielded).  (16) 

and  state  that,  in  effect,  inclusion  of  Tj  is  mathematically  equivalent  to  replacing  £iVU  by  EMf 

Tj- 

We  now  assign  the  "military  value"  V j  to  each  item  of  the  Jth  type.  (One  means  for  arriving 
at  the  military  value  is  addressed  by  BDM.I 11)  The  Vj  values  are  normalized  such  that 

2  y,=  1  •  (17) 

j=i 

The  Vj  is  a  linear  weighting  factor  to  the  Pdj  values.  Thus,  the  Pri  value  for  the  extended  target 
as  a  whole  is 

m 

P*=ZPvVi  •  (18) 

j^l 

where  shielding  has  been  considered  in  calculating  P,tJ.  Inserting  equation  (10)  into  equation  (18) 
and  simplifying,  we  find  that  Pd  is  given  by 


I  I 

P‘~2  +  2j2 


-  m 

?  'n  n 

-  [  j=i 


( T)E,i/ £r>0j )' 


Equation  (19)  can  be  compared  to  equation  (5).  the  probability  of  damage  to  a  point  target,  and 
we  obtain 

=  l  +  fn[(£r/ES0>'“]  .  (20) 

2  2v2 

"  BDM.  Inc.,  Theater  Nuclear  f  orce  Combat  C  apability  Degradation  Methodology.  Development,  and  Demonstra¬ 
tion.  BDM/W-78-  167-TR  (April  I97X).  (SECRET) 
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where  the  parameters  £50and  cr  represent  the  extended  target  by  a  single  entity  and  E,,  represents 
the  target/burst  characteristics  (weapon  yield,  distance  from  burst,  and  deployment  of  the  point 
targets  within  the  extended  target).  This  definition  of  £,.  means  that 

m 

(£„/£,„)''"  =  []  (T1ErifEwV*’>  .  (21) 

j=  I 

The  presence  of  Vj  in  equation  (21)  merely  serves  to  scale  the  value  of  the  corresponding 
<jj.  Thus,  in  addition  to  combining  the  transmission  factors  into  the  definition  of  £SUJ.  one 
could  similarly  incorporate  the  Vj  factor  into  the  definition  of  ct, ,(«-,.;  =  <jj/Vs).  That  is  to  say. 
the  vulnerability  parameters  may  be  easily  altered  to  include  the  effects  of  transmission  as  well 
as  military  value. 

To  calculate  £,/,  the  coordinates  and  yield  of  the  weapon  in  question  would  be  given,  as 
would  be  the  locations  of  the  sub-target  items  and  the  values  of  Vj.  r;r  Em.  and  Tj  for  each 
item.  Environment  algorithms,  developed  by  Sweeney  et  al.9  would  then  be  used  to  find  Erj  from 
equation  (II).  Then  the  value  of  the  right-hand  side  of  equation  (21)  would  be  calculated,  and  the 
result  inserted  in  equation  (20)  to  give  This  procedure  would  then  be  repeated  for  all  bursts, 
environments,  targets,  and  damage  categories.  This  methodology  provides  an  excellent  basis  for 
low-level  aggregation  of  targets. 

2.4  The  Vulnerability  Center 

Corresponding  to  Erj,  one  can  define  an  effective  distance  r,.j  that  can  be  determined  from 
the  E(r)  relations,  equation  (6).  The  distance  rrJ  is  not  exactly  given  by  the  distance  between  the 
target  geometric  centroid  and  the  burst  point,  but  depends  upon  the  particular  sub-target  distri¬ 
bution.  This  is  illustrated  by  considering  two  identical  point  targets  equidistant  (r,, )  from  a  burst, 
but  not  on  opposite  sides  of  the  burst.  Their  P„  values  are  the  same,  and  r,>  =  rri  is  located  on 
an  arc  of  radius  rrt  centered  at  the  burst  point.  Since  r ,,  =  r,,  is  located  at  the  intersection  of  the 
arc  with  the  perpendicular  bisector  of  the  two  targets,  r does  not  lie  at  the  target  centroid,  which 
is  not  on  that  arc.  Similarly,  irregular  target  shapes  are  cases  where  the  target  centroid  is  not 
colocated  with  r,..  The  interest  in  r,.  is  for  cases  when  its  location  can  be  closely  approximated 
by  the  target  centroid,  for  then  the  extended  target  can  be  replaced  by  a  point  target  located  at 
the  centroid.  To  calculate  the  effective  distance  rrj,  a  point-by-point  weighting  by  the  environment 
intensity  is  required.  Because  this  function  is  not  linear  in  distance.  rrJ  differs  from  the  target 
centroid  (which  is  found  by  a  linear  weighting  of  distances).  Because  we  are  concerned  with 
extended  targets  having  a  limited  extent,  it  may  prove  possible  to  approximate  rvj  by  the  target 
centroid.  We  now  consider  the  accuracy  of  such  an  approximation. 

If  all  items  are  of  the  same  sub-target  type.  Vj  =  1.  <tj  =  <r,  and  equation  (111  reduces  to 

n 

E„  =  I]  £/'*  .  ,22) 

i  =  1 
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where  the  subscript  j  has  been  deleted,  and  T  has  been  included  in  E-M.  For  all  environments. 
r,,  is  found  by  inserting  equation  (6)  into  equation  (22). 

n 

rfelr‘ •  =  [J  rf"‘elr'1'1  (23) 

i= i 

Thus,  r,,  is  independent  of  the  values  of  the  weapon  yield  and  the  vulnerability  parameters  of  the 
individual  sub-targets.  Because  V'and  Tjust  serve  to  scale  the  sub-target  vulnerability  parameters. 
r is  also  independent  of  V'and  T. 


From  equation  (23).  representative  of  a  single  sub-target  type,  the  value  of  r,  differs  from 

n 

that  of  the  arithmetic  mean  r  jn.  The  geometric  mean  r„  is  always  less  than  or  equal  to  the 

/  =  1 

arithmetic  mean  r„:M 


I  I  r! •  <  -  £  r, 


or 


r„  ^  r„ 


Then  the  bounds  on  r,.. 


II  r)‘*  r, 

"  "  ,r, 


or 


r  r  <  r 

'  ii  —  '  r  —  '  a  * 

may  be  demonstrated  by  a  two-step  process  based  on  equation  (23). 


(24) 


(25) 


(26) 


(27) 


We  rewrite  equation  (23)  in  the  following  form 


and  reduce  it  to 


(28) 


r*err-  =  r”t’Cr«  .  (29) 

In  the  first  step  in  the  process,  we  insert  r„  for  r,.  in  the  left-hand  side  of  equation  (29)  to  find 
that 


>•  B  /;  f  r y  7  ••  W  ^  f|| 

'  ft 1  =-'(/' 


(30) 


However,  since  r„  s  r„.  the  left-hand  side  is  smaller  than  the  right-hand  side,  implying  that  our 
choice  of  for  r,  is  too  small.  Thus. 

r„  s  r,  .  (311 


4  M  Abramowit/  and  I  A.  Stcgun.  editors.  Handbook  of  Mathematical  I  unctions  with  formulas.  Graphs,  and 
Mathematical  Tables.  National  Bureau  of  Standards.  Applied  Mathematics  Series  55  (June  |W»4). 
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In  the  second  step  in  the  process,  we  insert  r„  for  r,.  in  the  left-hand  side  of  equation  (29)  to  find 
that 

rf,e,  r-  l  />"•"  .  (32) 

Now.  since  r„  --  r„,  the  left-hand  side  is  too  large,  implying  that 

r,.^r„  .  (33) 

Combining  equations  (31)  and  (33),  we  obtain  the  desired  relation,  equation  (27). 


We  now  define  Ar,.  the  separation  of  the  individual  point  targets  from  the  target  centroid,  by 
the  scalar  relation 

r,  =  r,.  +  A  r,  .  (34) 

where  r,.  denotes  the  target  centroid.  If  A rjr,.  is  small  (about  10  percent  or  less),  then  r is 
located  at  the  geometric  centroid  (arithmetic  mean) 


to  within  an  error  that  is  given  by 


(35) 


£  (A r,/rr)2  .  (36) 

ZAI  /  =  1 

using  the  bounds  of  equation  (26).  The  negative  sign  of  equation  (36)  implies  that  the  vulnerability 
center  is  always  on  the  burst  side  of  the  centroid.  That  is.  the  closer  sub-target  items  are  the 
more  heavily  weighted.  With  a  target  extent  of  less  than  one-tenth  the  distance  front  the  target 
centroid  to  the  burst  point,  r,.  is  equal  to  the  centroid-to-burst  point  distance  to  within  a  fractional 
error  of  5  x  I  O'3.  This  small  error  mums  that  the  vulnerability  center  (the  location  of  the  point 
target  that  mathematically  replaces  the  extended  target)  is  well  approximated  by  the  target 
centroid. 


3.  APPLICATION  TO  BATTERY-SIZED  TARGET 

In  this  section  we  investigate  the  accuracy  i>f  the  approximations  developed  in  section  2  by 
applying  them  to  an  appropriately  sized  extended  target  (100  by  200  m)  that  is.  indeed,  composed 
of  discrete  points.  This  unit  is  the  nuclear-capable  155-mm  Field  Artillery  Firing  Battery.  Self 
Propelled.  The  prime  reasons  for  this  choice  are  that  an  idealized  deployment  (see  fig.  2)  of  this 
unit  has  been  developed,  the  sub-target  locations  have  been  identified,  and  this  unit  is  currently 
of  interest  in  applications  of  tactical  nuclear  warfare  analyses.1-  A  second  reason  for  the  choice 
of  this  unit  is  that  it  has  been  found12  that  this  battery  can  be  accurately  characterized  as  a  point 
target  located  at  its  "vulnerability  center"  when  the  nuclear  damage  mechanism  is  radiation- 
induced  incapacitation  of  the  unit's  personnel  (as  well  as  vehicle  overturn,  transient  radiation 
effects  to  radio  and  electronics  failure,  etc.). 

The  four  categories13  of  incapacitation  due  to  total  radiation  dose,  and  (he  vulnerability 


'*  ('.  K.  Spyropoulos  and  J.  Wicklund.  A  Method  lor  Assessing  the  Vulnerability  of  Small  l  nils  in  tactical  Nuclear 
Engagements  (II).  Harry  Diamond  Laboratories.  HOI .-TR- IKS I  (June  IWKi.  ((ONI  ) 

13  ('.  Stuart  Kelley.  Distribution  of  Nuclear  Probability  with  Distance.  Harry  Diamond  I  aboratoncs.  HDI  -TK-IXMt 
(August  IV7X). 
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Figure  2.  Deployment  of  the  field  artillery  battery.  155  mm  self  propelled. 


parameters  associated  with  these  categories,  are  shown  in  table  I.  The  criteria  for  the  incapaci¬ 
tation  of  the  personnel  will  vary  from  position  to  position  inside  the  unit's  boundary  because  the 
various  troops  perform  different  functions:  some  are  physically  demanding,  some  are  undemand¬ 
ing.  Table  2  shows  the  coordinates  (in  km)  of  the  troops  as  measured  from  the  centroid  of  the 
battery  [the  center  of  the  outline  of  the  battery  is  displaced  (-0.0075.  -0.0215)  from  the  battery 
centroid].  Incidentally,  for  this  example,  the  vulnerability  center  found  by  Spyropoulos  and 
Wicklund12  is  within  5  m  of  the  centroid. 

The  procedure  for  assessing  the  accuracy  of  the  approximations  of  section  2  is  to  detonate 
nuclear  weapons  at  various  distances  (R)  and  angles  (0)  from  the  battery  centroid,  perform  the 
exact  and  approximate  P,,  calculations  for  each  detonation,  and  compare  the  resulting  values  of 
P,, .  The  cases  considered  here  use  weapon  yields  of  W  =  1,  10.  and  100  kT  detonated  at  distances 
of  R  =  0.4.  0.5.  0.6,  0.7,  0.8,  1.0,  1.2.  and  1.5  km  from  the  battery  centroid  and  at  angles 
(measured  in  a  clockwise  sense  from  the  +  y  axis)  of  0  =  0.  45.  90.  135.  and  180  deg.  In  each  of 
these  cases,  we  consider  the  total  radiation  dose  experienced  by  each  of  the  21  groups  of  troops 
described  in  table  2. 


12  c.  h.  Spyropoulos  and  J.  Wicklund.  A  Method  for  Assessing  the  Vulnerability  of  Small  Units  in  Tactical  Nuclear 
Engagements  ( U ).  Harry  Diamond  l  aboratories.  HIM  .-TK- 1X5 1  (June  14710.  (CONI') 


The  total  radiation  dose  (in  rads  tissue)  experienced  by  the  ith  group  of  the  21  groups  is 

0,  =  8528  Wrf24"V“ 3  r,72r'  .  (37) 

where  r,  is  the  distance  between  a  detonation  and  the  /th  group  in  the  battery.  The  equation  for 
Pd  that  we  consider  to  be  exact  is  the  appropriate  form  of  equation  (I);  namely. 


The  number  of  troops  ///,  in  the  /'th  group  is  listed  in  table  2. 

We  first  consider  the  approximation  to  equation  (38)  that  occurs  when  erf(.i  )  is  replaced  by 
x.  Then  Pd  takes  the  form 

/>./.  =  ^|i'"^{1  'n(/>, //>,*»}  •  <39, 

Values  of  P and  P,n  were  calculated  for  combinations  of  W.  R.  and  H.  and  a  plot  was  made  of 
the  Plh  versus  P,ir  data  points.  Although  there  was  some  scatter  in  the  data  points,  they  are  very 
well  represented  by  the  solid  curve  of  figure  3.  The  dashed  line  on  that  figure  represents  P,„  = 
P(ll..  The  extent  to  which  the  solid  line  deviates  from  the  dashed  line  measures  the  extent  to 
which  P,„  deviates  from  P„,..  It  will  be  seen  from  figure  3  that  /\n  most  accurately  reproduces 
P,i,  near  PtU,  =  0.4,  but  diverges  widely  at  the  extremes  of  Pd,  =  0  and  1.0.  Over  the  range  of 
0.2  <  P,u.  <  0.8,  however.  Pdl  equals  P,t,  to  within  10  percent,  and  over  this  range  P,h  is  a 
useful  approximation. 

Consider  next  the  concept  of  the  effective  environment  intensity  ,  here  the  effective  dose.  £),.. 
to  the  battery.  The  effective  dose  is  defined  by  the  left-hand  side  of  equation  (21).  but  can  only 
be  determined  once  the  value  of  <r  has  been  found  for  the  area  target  as  a  whole.  The  quantity 
of  interest,  however,  is  the  right-hand  side  of  equation  (21):  applying  it  to  a  Pd  calculation  is 
exactly  the  same  as  calculating  P,,  by  equation  (39).  and  the  results  are  the  same  as  described  in 
the  previous  paragraph.  Therefore,  over  0.2  -5  P,,  <  0.8.  the  concept  of  the  effective  environment 
intensity  is  both  valid  and  accurate. 

The  remaining  portion  of  this  section  concerns  the  troops  in  category  2.  These  troops 
comprise  most  (63  percent)  of  the  battery's  troops:  therefore,  it  is  worthwhile  to  consider  the 
results  of  approximating  the  entire  battery  by  just  those  troops  in  category  2  alone.  As  a  first 
step  in  this  analysis,  the  exact  Pd  values  for  the  entire  battery  were  compared  to  the  corresponding 
exact  Pd  values  that  result  when  only  the  troops  of  category  2  are  included.  The  category  2.  Pd 
values  were  found  to  reproduce  the  exact  values  to  within  8  percent,  and  are  accurate  near  P d, 
=  0  and  1 .0. 

For  category  2.  the  effective  dose  is  given  by  the  appropriate  modification  of  equation  (22) 

D,  =  I]  (40) 

* -  I 


that  corresponds  to  an  effective  distance  r,  that  is  bounded  by  the  appropriate  modification  of 
equation  (26). 


II 


rm,i  5» 

•  i 


-  r<  -  75  2  m<r> 

'9  ,  =  i 


(41) 


For  the  cases  considered,  the  two  bounds  are  very  close  to  each  other  and  differ  by  less  than  I 


15 


/ 


Figure  3.  Accuracy  of  approximation,  equation  (2V).  to  exact  probability  of  damage: 
deviation  of  solid  line  from  dashed  line  measures  inaccuracy  of 
equation  <2S»)  in  reproducing  exact  result,  equation  (28). 


percent  of  the  value  of  r,..  Accordingly,  r,  is  accurately  specified  by  either  of  these  two  bounds. 
Also,  to  within  I  percent,  the  effective  distance  is  equal  to  the  separation  of  the  burst  from  the 
target  centroid. 

The  final  approximation  to  be  evaluated  consists  of  making  the  point  target  approximation 
using  D,  as  the  environment  intensity  ,  and  replacing  erf  v  by  x. 

( n(/>,  //)r,„2  >|  .  (42) 

The  results  of  the  calculations  are  quite  similar  to  those  shown  in  figure  3.  As  in  figure  3.  this 
approximation  gives  Pdl.  to  within  10  percent  over  0.2  <  s  0.8.  but  is  inaccurate  outside  this 
range. 


i  + 


\  2  <t  .2 


4.  SUMMARY 

This  report  deals  with  the  question  of  how  extended  targets  composed  of  discrete  points  can 
be  mathematically  treated  as  single  point  targets  when  the  analyst  calculates  target  damage  caused 
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by  a  nuclear  burst.  The  answer  involves  uncertainties  in  the  data  that  originally  gave  rise  to  the 
mathematically  complex  form  of  the  damage  function.  Given  these  uncertainties  and  the  damage 
function  values  that  are  of  interest,  the  damage  function  can  be  greatly  simplified,  thereby 
allowing  a  savings  in  the  calculations  needed  to  determine  damage  to  an  extended  target. 

The  procedure  for  calculating  damage  involves  sub-dividing  the  target  into  sub-targets  that 
have  similar  vulnerability  parameters.  These  parameters  can  be  easily  scaled  to  account  for 
transmission  and  "military  value"  factors.  The  resulting  formalism  permits  the  calculation  of 
"surviving  military  value"  of  any  such  target.  The  procedure  is  especially  appropriate  to  low- 
level  aggregation  schemes. 

The  accuracy  of  this  procedure  is  assessed  by  applying  it  to  the  problem  of  determining  the 
incapacitation  of  the  troops  of  a  155-mm  howitzer  battery  that  is  caused  by  a  nuclear  detonation. 
Consistent  with  the  results  of  other  analyses,  the  concept  of  effective  distance  is  shown  to  be 
appropriate,  and  therefore  this  extended  target  can  be  considered  as  a  point  target  that  is  located 
at  the  centroid  of  the  extended  target.  The  simplified  damage  function  is  found  to  have  acceptable 
accuracy  for  most  damage  probabilities. 
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ATTN 

CHAIRMAN,  EDITORIAL  COMMITTEE 

ATTN 

SOKOLOSKI,  M.  M. ,  00210 

ATTN 

PEPERONE,  S.  J.,  36000 

ATTN 

SANN,  ] 

K.  H.,  11100 

ATTN 

CHIEF, 

20240 

ATTN 

CHIEF, 

21000 

ATTN 

CHIEF, 

21100 

ATTN 

CHIEF, 

21200 

ATTN 

CHIEF, 

21300 

ATTN 

CHIEF, 

21400 

ATTN 

CHIEF, 

21500 

ATTN 

CHIEF, 

22000 

ATTN 

CHIEF, 

22100 

ATTN 

CHIEF, 

22300 

ATTN 

CHIEF, 

22800 

ATTN 

CHIEF, 

22900 

A.TN 

ISLER, 

W.  E.,  13500 

ATTN 

SHARE, 

S.,  22800 

ATTN 

BELLIVEAU ,  L. ,  22100 

ATTN 

WICKLUND,  J.,  22100 

ATTN 

SPYROPOULOS,  C. ,  22100 

ATTN 

FEMENIAS,  R. ,  22100 
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